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CN , Abstract 

o 

' We evaluate Polyakov loops and string tension in two-dimensional QED with both massless 

and massive TV-flavor fermions at zero and finite temperature. External charges, or external 
| electric fields, induce phases in fermion masses and shift the value of the vacuum angle parameter 

9, which in turn alters the chiral condensate. In particular, in the presence of two sources of 
opposite charges, q and —q, the shift in 9 is 2n(q/e) independent of N. The string tension has 
O 1 a cusp singularity at 9 = ±n for N > 2 and is proportional to m 2Ar /( Ar + 1 ) at T = 0. 

<D " 

> 

^ . Two-dimensional QED, the Schwinger model, with massive iV-flavor fermions resembles four- 

dimensional QCD in various aspects, including confinement, chiral condensates, and 9 vacua Q- 11]. 
Much progress has been made recently in evaluating chiral condensates and string tension in the 
massive theory |12|-[p!q]. In this paper we shall show that the three phenomena, confinement, chiral 
condensates, and 9 vacua, are intimately related to each other. In particular, the string tension in 
the confining potential is determined by the 9 dependence of chiral condensates {ipip). 

The behavior of the model is distinctively different, depending on whether N = 1 (one-flavor) 
or N > 2 (multi-flavor), and on whether fermions are massless or massive. The massless (m = 0) 
theory is exactly solvable. ( Tptp ) e / for N = 1, but ( i> i/j ) e = for N > 2. 0, |§ In either 
cases the string tension between two external sources of opposite charge vanishes f|, [T^] . 

In the massive (m ^ 0) theory (ijj'ip)g is proportional to m( N ~ l "( N+ ^ cos 2N >( N+1 '(9/N) at 
T = where 9 = 9 - 2vr[(6» + vr)/2vr] §, |||. For N > 2 the dependence on m is non-analytic. It 



also has a cusp singularity at 9 = ±7T. A perturbation theory in fermion masses is not valid at low 
temperature. 



The confinement phenomenon can be explored in various ways. One way is to evaluate the 
Polyakov loop P q {x) = exp {iq j£ dr A (t,x)} at finite temperature T = /T 1 @-p|. F(T) = 



—Tin (P q (x) ) or —Tin { P q (x)^ P q (y) ) measures the increase in free energy in the presence of an 
external charge q or a pair of charges q and —q. In particular, the latter is written as a\x — y\ for 
large \x — y\ where a is identified with the string tension. This method has the advantage of giving 
the temperature dependence directly. 

Alternatively one may determine the ground state to evaluate the change in the energy density 



(at T = 0) when a pair of sources of charge q and — q is placed ||, g4|, 25 1. This method has the 
advantage of showing how external charges affect the 9 parameter and chiral condensates. 

We employ both methods in a unified manner. Years ago Coleman, Jackiw, and Susskind 
showed the confinement of fractional charges in the N = 1 theory adopting the latter method 
01. Recently Hansson, Nielsen and Zahed applied functional integration method to evaluate the 



Polyakov loop correlation function [12]. The argument has been generalized to finite temperature 



by Grignani et al. [15]. Ellis et al. [|| and Gross et al. 14] have presented the mechanism of 
confinement in terms of soliton solutions in the bosonized form. All of these arguments are given 
in the one flavor (N = 1) case and rely on the validity of a perturbation theory in a fermion mass. 

Recently chiral condensates with arbitrary fermion masses m, vacuum angle 9 and temperature 
T have been evaluated in the A-flavor model [jDJ. The problem was reduced to solving a quantum 
mechanical system of A degrees of freedom. It was shown that the m — > and T — > limits do not 
commute for A > 2. In particular, the m-dependence of physical quantities is singular at T = 0. 

We analyse the model 

N 

_ 77 WW _i_ 



mass 

a=l 

N 

£mass = " E m «{ ^ M - + e ~ iSa M a } (™« > °) 

a=l 

M a = i> Ul-^a (1) 



defined on a circle with a circumference L [24|-|35]. The model defined at finite temperature 



[46], on a torus or sphere [47|-|52], or on a lattice or light-cone 53 ]-[ 52] , has been also extensively 
discussed in the literature. 

We impose boundary conditions A^(t,x + L) = A^(t,x) and ip a (t,x + L) = —ip a (t,x). On a 
circle the only physical degree of freedom associated with gauge fields is the Wilson line phase 

e w (t): 



3 »6w(t) = exp , ie 



jie J dxAi(t,x)^ . (2) 



In Matusbara's formalism the finite temperature field theory is defined by imposing perodic or 
anti-periodic boundary conditions in the imaginary time (r) axis on bosons or fermions, respectively. 
Mathematically the model at finite temperature T = /3 _1 is obtained from the model defined on 
a circle by Wick rotation and replacement L — > (5, it — > x and x — > r. Furthermore the Polyakov 
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loop of a charge q in the finite temperature theory corresponds to the Wilson line phase: 
Pq(%) = exp liq J dr Aq(t, x)| <==^- exp ji - ©w(0 



(3) 

We bosonize fermions in the Coulomb gauge || |l7], 24, 13]. Take 7 M = ((Ti,io"2) and write 
ipj = (^2. , ) . In the interaction picture defined by free massless fermions 



}_ qo, e ±i{q%+2-Kp\(t±x)/L} . e ±iy/A^<f, a ± {t,x) . 



e 2mp ± |phys) = |phys) 
where C% = exp{ivr (p b + + pt)} and = exp{^ ELi(p+ ~ P-)}- Here 

[4>p±] = 

oo 

= ^(^-^{c^^e-^^ + h.c.} , [4, n ,4!m]=^nn 

71=1 



(4) 



(5) 



The : : in (||) indicates normal ordering with respect to (c n ,cjj. In physical states p± takes an 
integer eigenvalue. 

Conjugate pairs are {p a ,q a } = + P-),q+ + q"}, {PaAa} = {p+ - j£,§(«+ - ql)}, 

{Pw, @w} ; and {n a , a = 0" + (j}°L\. The Hamiltonian in the Schrodinger picture becomes 

H tot = H + Htf, + H mass + (constant) 

e 2 L 2 iV r A 2vrJ^ ^ 2vr r ^ 2tt A r 2 1_ 2 -, 

rL 1 E{n^ + (tf} + 7(E4 



a=l 



a=l 



a=l 



^mass = f L dxY,m a { e i5a M a + e -«» M] } . 

•/ „ 



In the mass term 



i\/4:1T<j>a] 



(6) 



(7) 



where N^[- ■ •] indicates that the operator inside [ ] is normal-ordered with respect to a mass /x. In 
general a mass-eigenstate field Xa with a mass is related to <p a by an orthogonal transformation 
Xa = U aa (j)a- In (§) we have 



No[e' 



N 



B a =l[ B{ii a L) 



a=l 

(U aa ) 2 



B{z) 



a=l 

Z 

4tt 



exp <^ 7 + 



7T 



du 



z h (e uz - l)Vu 2 - 1 



(8) 
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As \p a ,H tot ] = 0, we may restrict ourselves to states with p a = 0. 

Hmass gives rather complicated coupling between the zero and (p a (xa) modes, whose effects are 
non-perturbative for N > 2. As in previous papers [13| the vacuum wave function is written in the 
form 



where 



r2n N-i 

$vac(#) ) = / ]Jdip a f((p a ; 9 cS ) |$ (^a + S a ; 6) \ 

J ° a=l 



|^0(^;6»)> = (27T)~ Ar / 2 e in0+i ^a=i^ ^(n+r 1 ,-,n+r N - 1 ,n) ) 

{n,r a } 

9 N 

<ew,?a|^ Bl, - ,Bw) ) = «o(9w + f E"-) {2^)- N/2 e^ a ^ aqa 

a=l 

u 



(9) 



N 

e eS = e-Y,Sa . (io) 

a=l 

We have generalized the expression to incorporate the phases 5 a 's in the mass parameter in (|]). 
The ansatz above for the vacuum is good for m a <C e. 

The eigenvalue equation (Hq + H mass )\Q VSiC (6) ) = E\& vaiC (9) ) reads 



where 



{ - + V N {<p; 9 cS )} f{<p; 6 cS ) = e /(</>; 8 cS ) (11) 
d 2 2 d 2 



A 



NL N 
V N (^0 eS ) = - e -«/W Vm a B a cos^ 

(N - 1)tt ^ 

N-l 
a=l 

and e = NEL/2ir(N - 1). 

For the x« fields, the vacuum is defined with respect to their physical mass (jL a 's which needs 
to be determined self-consistently from the wave function in (|9|). In the symmetric case m a = m 
one has /i 2 = • • • = /xjv, Ml = ^ + M2 and ^ a = B(^iL) l / N B^L)^^^ = B. The potential is 
reduced to 

JV 

Vn(<P, Oes) = ~«Q X! COS( ^a 
a=l 

_ iVmL /JVmL 

Further /i| is determined by 

Ml = ^=^e-^(cos^) / (14) 
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where the /-average is given by (g((f) ) f = J [dip] g(ip)\f(ip)\ 2 . We have made use of the fact that 
( e lifa ) j is independent of a. Eqs. (|ll|) and ( |i~4| ) are solved simultaneously. Evaluation of these 
equations was given in |13[ |. An important point in the following discussion is that /(</?; # e ff) or 
^vac(^) is determined solely by m, L (or T), and e ff- 

Now let us evaluate the Polyakov loop at finite T, or equivalently ( e * fc0 w(i) ) e Q n a circle. The 
parameter corresponds to q/e where q is the charge of an external source. Since the expecta- 
tion value is i-translation invariant, it is sufficient to evaluate at t = 0. Making use of @, one 
immediately finds 

( <&(0'; ^')|e lfcew <p) ) = <? 27r (0 - 9> - 2-nk) Y[ 6 2 „(<p a - <p' a - ^) • e ~^ 2 ^/4N (15) 



a=l 



It follows that 

(Pke)e,T = (j keW )9,L=T-i 

' for k 7^ an integer 

e~ k2 ^/ iNT J[dip]f(<p a ;9 eS yf(<p a + ^;e eS ) for fc = an integer . ( 16 ) 

The vanishing of the Polyakov loop for a fractional A; is due to the invariance under large gauge 
transformations. The Hamiltonian (||) is invariant under 0\y — ► ©w + 27T and p a — > p a — 1. In 
other words, 

A 7 " 

{/ = exp (2niPw +i^q a ) 
[U,H tot ]=0 

U\$ V3C (e)) = e ie |$ vac (#)> , (17) 

which implies the vanishing of ( e lk&w )g for a non-integer k. 

In the JV = 1 (one-flavor) case, there is no (p a degrees of freedom. Eq. (|l~6| ) reduces to ( Pk e )bt = 
e -k 2 nfi/4T ^ w j 1 j c j 1 agrees with the result of Grignani et al. The factor k 2 irfi/4 is understood as the 
self-energy of the source [[0| . 

In the N >2 (multi-flavor) case the overlap integral for the f(<p;0 e s) factor becomes relevant. 
In the massless (m = 0) case, however, f(ip) is constant as the potential V^{ip; 6 e s) in ( |TT| ) vanishes. 
Hence ( P^ e )qt = e~ k2nfl ^ 4:NT in the massless theory. 

When m ^ 0, the overlap integral needs to be evaluated numerically. In two limits, namely 

N 

T//J, <C (m/fi) N + 1 and T//i> 1, analytic expressions are obtained. It is instructive to examine the 
free energy F e (T) = -Tln(P e ) e>T for k = 1. At sufficiently low T = L~ l , k > 1 in (||) so that 
f(<p) has a sharp peak at the minimum of the potential 



VT l = -S > 0* = OcE - 27T • (18) 



6» cff + 7T 



2tt 



/(</?) is approximately given, up to a normalization constant, by 



f = expi - \ l K COS — I }^ t fa+ }^ VaPb J f ' <Pa = <Pa ~ ■ (19) 

I ' ^ a=l a<£> ' J 
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N=3 m/n =0.01 




+ Overlap 


/ < P > 


integral 

+ 




+ o 


- 

(a) 


1 minirt ft-akjjLg a a a ih a, ft i^i . i 1 1 l. 





0.1 



F/ji 



N=3 m/n =0.01 



<>P x I / >/H 



(b) 



T/ > 

Figure 1: (a) The Polyakov loop (P) and the overlap integral in (16) are plotted as functions of 
Till for k = 1 (q = e), N = 3, mj \i = 0.01 and 6 = 0. (b) The free energy F/fi = -(T/fi) In ( P) 

and the chiral condensate (ipip) I M are plotted with the same parameter values. 



Hence the overlap integral gives an additional damping factor in ([16]). In the opposite limit T//j> 
1, kq <C 1 so that / ~ constant. Hence we find, for an integer k = q/e, 



F q {T) 



2 ( ~ JV 

k J i i (at i\f-> i m 6leff ^ n+i 

_|l+(iV-l)^-co S — 

A; 2 7r/x 
4A^ 



for T < m N + 1 /j, N + 1 
for T 3> // . 



(20) 



The free energy is finite. It does not diverge even at T = 0. An integer external charge is screened. 

In fig. 1 we have depicted the temperature dependence of the Polyakov loop, free energy, and 
chiral condensate in the N = 3 case with m/fj, = 0.01 and 9 e Q = 0. Notice that all these quantities 
show cross-over transitions, but at different temperature. 

The vanishing of the Polyakov loop for a fractional charge q = ke does not necessarily implies 
the confinement as it follows from the gauge invariance. To obtain information on the confinement 
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or the string tension, we evaluate the Polyakov loop correlator: 

G q (x) = ( P q (x)^P q (0) W ^ G q=ke {t) = < T[ e - jfc0 wW e **ew(0) ] ) L . (2 1) 

Without loss of generality we suppose that x > and t > (x *-+ it). Note that G q {t) is gauge 
invariant. 

In the N = 1 model the correlator G e (t) was first evaluated by Hetrick and Hosotani For 
a general (/, G q {x) has been evaluated by Hansson et al. [12] and by Grignani et al. [15] . Consider 

G(t; k,l) = ( T[ e ~ ifce w(<) e +^e w (o) j } & L ( 22 ) 

in the massless iV-flavor model, which we denote by G^°\t;k,l). In this case the zero modes 
(6w, ?o) and oscillatory modes decouple so that the model is exactly solvable. The Heisenberg 
operator 9w(t) can be expressed in terms of Schrodinger operators 

8w(*) = 0\ycos/zi + — jy-Pw sin fit + — ^2p a (cos fit - 1) . (23) 

a=l 

Making use of (^), we find 

e - lkew (t) e +«Gw(o) = exp | _ sinMt + !^ sin ^ {k cos ^ t _ /) 

x exp | ^— (cos fit — 1) ^2 Pa — i(k cos fit — i)B\v| 

f ikir/iL 1 , 

x exp < — — sin^tPw^ . (24) 

In taking the vacuum expectation value of (^), we encounter the factor 62^ [0 — 0' + litik — I)] as 
in (|l5|). It follows that 

for k — I ^ an integer 

G(°)(t;A;,0 = -i exp |_ ^ (fc 2 + p _ 2We -W[)J for - Z = an integer. ( 25 ) 



This implies that the increase in the free energy in the presence of a pair of charges q and —(7 is 

F pair (r) (0) = _ r i n (p 9 ( x )tp g(0))e = ^0 2 (1 _ e -MN) . (26) 

In the massless theory external charges are completely shielded and the string tension vanishes as 
shown by various authors |, ||, 

If fermions are massive, the situation qualitatively changes. In the literature only the N = 1 
case has been analysed for which a perturbation theory in fermion masses is valid. Restricting 
ourselves to N = 1 with k = I and 5 = (therefore 6 e g = 9), we find 

G^ (t; fc, fc) = G(i; fc, fc) - (t; fc, fc) 

= -i /°° (JsjlTle-^O'e^O'ffats)']) - G( >(t; fc, A;) < #int(s) J >} 
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Hint{s y = _ll^^^ L dx { e ^)'^ . (27) 

Here the superscript / indicates the interaction picture defined by a massless fermion. To O(m) 
the (j) field part of Hi nt does not contribute. 

In the second term in the expression of G^\ {Hi nt (s)) = mBe~ n ^ L cos 8. In evaluating the 
first term we need, in addition to (|23|), 



x I dsle 



T 2 

q(t) =q+—jr (G w + 27rp)sin^ + 2vrP w (l -cos/ut) . (28) 

A useful identity is 

{e ± ig ( s y e -ike w (tV e *e w (o)' )(j = G V)(t;k,k) ■ e* ie e~^ L e ^k(i-e^) e -^ _ (29) 

Without loss of generality we take t > 0. The integral over s in (^7|) splits into three parts: 
/°oo> Jo> anc ^ /t°°- ^ ^ s eas y to cri eck that the first integral is the same as the third integral after a 
change of variables, and each of them vanishes. The manipulation is justified with the hypothesis 
of adiabatic switching of interactions implicit in the derivation of Gell-Mann-Low relations. 

The second integral gives the sole contribution to G^: 

G (1) (i; k, k) = imBe- n/f " L G ( - 0) (t; k, k) 

>^k e - t nk(e^(^ + e-^) _ ^ + {& ^ _^ fc ^ 1 (3Q) 

The integral is expressed in terms of Bessel functions. The correction to the free energy is, after 
making a Wick rotation it = x > 0, 

FP air (x,T)W = -TG W /G i0) 

= -2|x|mr5e _7r/Mi |j (27rA;z)cos(6'-27rA;) - cos 6 

+ —- Y e~ i( - 6 - 27Tk) + (-\y e +W~2vk) -n _ z n\j( 2 Trkz) } (31) 

W\ V / n J 

where 2 = e A*! 33 !/ 2 . For > 1, z C 1 so that 

F 9 pair (x,r)W ~o\x\ 

a N=1 = -2mTB(^)e^ T ^{co S (9-^)- C o S e} . (32) 
Here a is a "string tension". Since {1^1^)0 = — 2Te~ TT ^B(fj,/T) cos 8, we find 

=1 = m{(lp1p) e _ 2 ^ g/e) - <^>e} • (33) 



In other words, the major effect of a pair of external sources of charges q and — q is to shift the 9 
parameter in the region bounded by the sources by an amount 2n(q/e), which changes the chiral 
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condensate [T^] . A linear potential results because of this. We shall show below that this is true 
even for N > 2. The expression ( |33] ) is valid at arbitrary temperature. 

The string tension a can be either positive or negative, depending on the values of and q/e. 
This implies that the 9^0 vacuum is unstable against pair creation of sufficiently small fractional 
charges. 

The perturbation theory in fermion masses cannot be employed in the N > 2 case as physical 
quantities are not analytic in m at T = [||, The perturbation theory can be applied only in 
the high temperature regime. 

There is a better way to explore the problem. We place external charges on a circle and solve 
the Hamiltonian as was done in the N = 1 case in and pq ]. 

In the presence of external charges C ext = — ^oPext- Gauss's law implies 

d x E ext (x) = /wO) • (34) 

Let us restrict ourselves to static sources p C xt(x) where J$ dx p cxt = 0. Then E ext (x) = E^ xt — 
A^ ± {x) 1 . Here E^ t is constant and A§ xt (x) = — Jq dy G(x — y)p ey ±{y). In particular, for a pair of 
sources located at x = and at x = d. 



Pey±{x) = q{5 L (x) - 5 L (x - d)\ 



E cxt (x) = E^ t + E cxt (x)W 



(l — 4) for < x < <i 



= -Ar(x)'=\ d L ' (35) 
-q— for d < x < L 

1j 



Note that J L dx E ext (x) w = 0. 

Suppose that m a = m <C p. The total charge density is jg M = J2a=i ^t^a + Pcxt, and the 
Coulomb energy becomes 



1 f L 

^Coulomb = ~o / ^^EM^M^ ~ ^EM^) 

2 Jo 

= f L dxUpxi-E^f . (36) 
JO 2 

Here xi = N^ 1 / 2 J2a=i <fia- In view of (|36|) we write 



Xi = X\ + Xl 
d 2 



dx 2 

The total Hamiltonian is now 



^tot W — Hq + H c [ + H x + H-a 



(-i + M 2 )xf = M^ • (37) 
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H x = £ dx\ : jll? + (xl) 2 + f X \ + E(n 2 + x'J)} : 



(38) 



Hq and H mass are given by (|). When m < /i, jUi ~ //. In 



In other words, the net effect of a pair of external sources is to give x-dependent fermion mass 
phases 5 a = ^Air/Nxi{x) in (||). Suppose that yT 1 <C d <C L. Sufficiently away from the source 

^l(l-i)=^ for -i <x<d _ fl -i 

_fIlA= S ont for d+ -1 <x<L -^-l . 

ive L 



Finding the exact form of the ground state wave function of (38) is rather involved. Instead, we 
content ourselves with finding an approximate wave function, noticing that <5^ ff is almost constant 
between the two sources. 

The entire circle is divided into the two regions, the inside region < x < d and outside region 
d < x < L (/U _1 <C d <C L). For the evaluation of local physical quantities in each region, one 
can approximately write the ground state as a direct product of ground states in the two regions: 
\^g) ~ l^)in ® l^out- I n the absence of sources, |\P S ) ~ |0) in ® |0) out . In the presence of sources 

\* g )~\8 + 59;5 i ») by ®\9 + 59;5™ t ) out . (41) 

In addition to the effect of <5^ ff an overall shift 59 in the 9 value results as <5° ff is x-dependent. After 
all there is only one 9 parameter globally. 

To determine 59, we utilize the fact that local physical quantities in the infinite volume limit 
L — > oo must reproduce results in the Minkowski spacetime. In particular, physics in the outside 
region d + fi^ 1 < x < L — fi^ 1 (L — ► oo) must be essentially the same as physics in the absence of 
sources. In other words \8 + 59;5 ont ) out ~ |#;<5 out = 0) out . As shown above physics depends on 6 
through the combination 6 c r = — J2a=i This determines 

59 = N5 out = • - . (42) 

e L 

Hence, the net effect is summarized by 

\*g) ~ |^eff)i„ ® \0) out 

e eS = e-^- ■ (43) 

e 
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Consequently, the change in the energy due to the external sources is, to 0(d/L), 

AE = Nmd{(^i;) ecS -{^i ) ) e } (44) 

so that the string tension is 

a = Nm{(^tlj) 9 _ 2n{q/e) - (ipip) e } • (45) 

The result generalizes to finite temperature, as was seen above in the N = 1 case. 

Note that the parameter 9 is not completely equivalent to the electric field E. Indeed, in 
the absence of sources (E) e = {eP w -A' ) e = 0. External sources, or external electric fields, 
induce effective fermion mass phases 5^ in the Hamiltonian, which in turn changes the effective 
9 e g through the chiral anomaly. We also remark that the Coulomb energy ( |36| ) is 0[(d/ L) ]. The 
linear potential results from the change in the chiral condensate. 

The chiral condensate (tp , tp)g at arbitrary temperatue T = L^ 1 has been evaluated in ref. [|b|]. 

With given m, the dependence of a on charge q is essentially the 9 dependence of ( ip if) ) g . At T = 
it has a cusp at # c fr = tt (mod 2tt). More explicitly 

at / x 2AT 

T=Q 2 ^ I ~ ~ 171 » N+1 




Notice the singular dependence of a on m as well. A mass perturbation theory cannot be employed 
at low temperature for N > 2. 



In the high temperature limit 



tt(JV-I) 




,_ . 2N for m N + 1 fi N +i CT«/i / x 

{i/ji/j) 9 = — — -m { p p (47) 



for T>/j 



for N > 3. There appears no 9 dependence to this order. Hence the string tension a is at most 
0(m 3 ) in this regime. 

For N = 2, the expressions for (^V)e m ( P^D * s multiplied by a factor 2 cos 2 \9. Therefore 

1 , I ^— for mkuh <C T < u 



7T 



e ~7rT/ M for T > /i 



For a general value of T/jU, a must be evaluated numerically. In fig. 2 we have displayed q/e 
dependence of a / > 2 at T/fi = 0.003, 0.01, 0.03, and 0.1 with m/n = 0.01 and 9 = in the iV = 3 
case. One can see how a cusp behavior develops at q/e = 0.5 (9 c s = it) as the temperature goes 
down. At higher temperature the magnitude of the string tension rapidly diminishes. For instance, 
aji? = 1.8 x 10~ 7 at T/fi = 1 and q/e = 0.5 . 

In this paper we have shown that the confinement of fractional charges in the massive N- 
flavor Schwinger model results from the effective change in the 9 parameter which alters chiral 



11 



N=3 m/ji =0.01 








/x 


\ 


T / (j, 
T/ji 


= 0.003 - 
= 0.010 

- 


//''T/u 


= 0.03 








= 0.10 







q/e 



Figure 2: String tension a/// 2 in (45) is depicted as a function of g/e at various temperature in the 
N = 3 case with mj\i = 0.01 and = . A cusp develops at q/e = 0.5 in the T — > limit. 



condensates. In the multi-flavor case (N > 2) the string tension at zero temperature has singular 
dependence on fermion masses and the ^-parameter. 
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